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Abstract. Following a first paper on this subject (Ciotti 
1991, hereafter Paper I), we study the dynamical prop- 
erties of spherical galaxies with surface luminosity profile 
described by the i? 1/,m -law, in which a variable degree of 
orbital anisotropy is allowed. The parameter m for the 
present models covers the range [1, 10]. For these models 
we study the self-consistently generated phase-space dis- 
tribution function (DF), and we derive - as a function of 
m - the minimum value of the anisotropy radius for the 
I model consistency (i.e., in order to have a nowhere nega- 
■ tive DF) . Then we study the region in the parameter space 
\ where the i? 1 /™ models are likely to be stable against 
radial-orbit instability, and we compare its size with that 
of the larger region corresponding to the consistency re- 
quirement. For stable anisotropic models the spatial and 
projected velocity dispersion profiles are obtained solving 
the Jeans equation, and compared to those of the glob- 
ally isotropic case, already discussed in Paper I. The rel- 
evance of the results in connection with the Fundamental 
Plane (FP) of elliptical galaxies is pointed out: the effect 
on the projected velocity dispersion due to the maximum 
orbital anisotropy allowed by the stability requirement is 
well within the FP thickness, and so no fine-tuning for 
anisotropy is required. Finally, the Velocity Profiles are 
constructed as function of the projected radius and for 
various degrees of anisotropy, and their deviations from a 
gaussian discussed. 

Key words: galaxies: elliptical - galaxies: kinematics and 
dynamics - galaxies: structure of 



1. Introduction 

The i? 1 / 4 -law [Eq. (1) below, with m = 4] was introduced 
by de Vaucouleurs (1948) to describe the projected lu- 
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minosity density (or surface brightness) I{R) of ellipti- 
cal galaxies, and has worked remarkably well. It has no 
free parameters and depends on two well defined physi- 
cal scales: a characteristic linear scale, i? e , and a surface 
brightness factor, I . 

A natural generalization of this empirical law was first 
proposed by Sersic (1968), as the i? 1 / m -law. From an ob- 
servational point of view, the i? 1 / m -law has been widely 
used (see, e.g., Davies et al. 1988; Capaccioli 1989; Makino 
et al. 1990; Young & Currie 1994; Andredakis, Peletier & 
Balcells 1995; Courteau, de Jong & Broeils 1996). In par- 
ticular, for normal ellipticals and brightest cluster galax- 
ies, a correlation between their luminosity L and the value 
of m has been found, in the direction of m increasing with 
increasing L (Caon, Capaccioli & D'Onofrio 1993; Graham 
et al. 1996). 

From a theoretical point of view much less work has 
been done on the i? 1//m -law, its apparent universality, and 
its possible applications to the problem of the FP of ellipti- 
cal galaxies (Djorgovski & Davis 1987; Dressier et al. 1987; 
Bender, Burstein & Faber 1992) and only one- component, 
spherically symmetric, globally isotropic models have been 
studied (Paper I; Hjorth & Madsen 1991; Ciotti, Lanzoni 
& Renzini 1996; Graham & Colless 1996). Considering the 
extensive use of the i? 1 / m -law we plan to extend the in- 
vestigation of this class of models. 

The interest in the study of the dynamical properties of 
the R 1 ' m models is also renewed following recent ground 
based observations (M0ller, Stiavelli, & Zeilinger 1995), 
and Hubble Space Telescope ones, showing that the spatial 
luminosity distributions of elliptical galaxies approach the 
power-law form p(r) cx r -7 at small radii, with < 7 < 
2.5 (Crane et al. 1993; Jaffe et al. 1994; Ferrarese et al. 
1994; Lauer et al. 1995; Kormendy et al. 1995; Byun et 
al. 1996; de Zeeuw & Carollo 1996). Since the deprojected 
density of the R}l m models increases toward the center as 
r -(m-i)/m f or m > 1 (Paper I), this family of models can 
be used to study power-law galaxies with < 7 < 1. 
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equals half of the total luminosity). The defining parame- 
ter is to, a positive real number, and 6 is a dimcnsionlcss 
parameter whose value is determined by the definition of 
R c . The function b = b(m) [Eq. (P5)] is very well fit- 
ted by the linear interpolation b(m) = 2m — 0.324, for 
0.5 < m < 10, with relative errors smaller than 0.001, and 
6(4) = 7.66924944 (Paper I). The total luminosity L m is 
given by L m = I R 2 C x L m , with 



2nm 



r(2m), 



In particular, in this paper we study one-component, 
spherically symmetric, anisotropic i? 1 /" 1 models, in which 
orbital anisotropy follows the widely used OM parameter- 
ization (Osipkov 1979, Merritt 1985). We numerically con- 
struct their phase-space DF, and we determine the region 
in the parameter space where anisotropic R 1 '" 1 models 
are consistent, i.e., their DF is positive over all the acces- 
sible phase-space. We then investigate the models stability 
against radial orbit instability, by using the global stabil- 
ity parameter, comparing the radial and tangential kinetic 
energies. In this way we approximately bound the region 
in the parameter space where the anisotropic models are 
consistent but (with high probability) unstable. 

Having determined the region in the parameter space 
where the R 1 /" 1 models are physical, we study their main 
properties. First of all the spatial and projected velocity 
dispersions are derived and discussed. Then we construct 
their Velocity Profiles (hereafter VPs), and fit them using 
the Gauss-Hermite series, discussing their deviations from 
a pure gaussian for different to and anisotropies. 

In a third paper (Paper III, Ciotti & Lanzoni 1996) 
we study in detail the properties of the DF, VPs, and 
velocity dispersion profiles of two-component spherically 
symmetric R x l m models, in order to understand how the 
superposition of a dark matter halo modifies the internal 
dynamics of the models and their observational properties. 

The paper is organized as follows. In Sect. 2 the basic 
properties of the R}/ m models, already discussed in Pa- 
per I, are summarized. In Sect. 3 the DFs for anisotropic 
models are derived and discussed, together with the limits 
imposed on orbital anisotropy for the models consistency 
and stability. In Sect. 4 we obtain the velocity disper- 
sion profiles of some characteristic models, and some ob- 
servationally interesting properties are presented. A short 
discussion on the implications of the obtained results on 
the problem of the FP of ellipticals is given. In Sect. 5 
the models VPs are constructed, and then fitted using 
the Gauss-Hermite series, thus quantifying the departures 
from gaussianity. Finally in Sect. 6 the main conclusions 
are summarized. 

^ (s) ^^[V 2 ;(";-W exp (- 6s v TO)s a 

2. A summary of the properties of spherical R 1 ' 171 2mb m 
models 



(2) 



where T is the complete gamma function [Eq. (P4) , Erdely, 
Magnus, Oberhettinger & Tricomi 1953, hereafter EMOT, 
vol. I, p. 1]. 

The most important deprojected quantity associated 
to I(R) is the luminosity density v, which is related to 
the mass density via p(r) = Tv{r), where T is the mass- 
to-light ratio, and r is the spatial radius. We assume a 
constant T, so that the main quantities (mass inside r, 
potential, velocity dispersion, etc.) depend only on the 
luminosity density v(r), which is related to the surface 
brightness profile by an Abel integral equation (see, e.g., 
Binney & Tremaine 1987, hereafter BT): 



v(r) = -- 



f 

J r 



dl dR 



dR^W 



= — lim 



I(R) 



Vr~ 2 



(3) 



The second term in the r.h.s. of Eq. (3) is zero for any 
positive value of to, and the resulting luminosity density 
v mif) = (Io/R c ) x v m (s), where s = r/R c , is extensively 
discussed in Paper I. The asymptotic behaviour of v m for 
r — > oo is given in Eq. (P8), while for r — ► one obtains: 



v m (0) = — r(l-m), 

7T 



v m {s) ~ - In 

7T \0S 



TO < 1 



TO 
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(4) 



(5) 
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Here we give a short summary of the basic properties of 
the R}l m models, as derived in Paper I. A full treatment 
is given there, and we will refer to equations in that paper 
as Eq. (P#). 

The spherical R 1 ^" 1 models are defined as a one- 
parameter family of stationary stellar systems, with sur- 
face brightness profile given by 



I(R) = I cxp(— brj 



l/r, 



(1) 



for to > 1, and where B(x,y) is the complete beta func- 
tion [Eqs. (P9)-(PA5), EMOT, vol. I, p. 9]. It should be 
noted that for to > 1 the density diverges at the origin as 
r 0—m)/m. therefore the divergence is worse for higher-m 
models. Finally, as in Paper J, we consider the relative po- 
tential ip m (r) — GTI R e x ip m (s), where G is the gravita- 
tional constant. Unfortunately ip m {s) cannot be expressed 
in terms of elementary functions, but, at variance with the 
density, for r = it converges for all to: 



where I is the central surface brightness, 77 = R/R c , R 
is the projected radius, and R e is the effective radius (i.e. 
the projected radius inside which the projected luminosity 



V>m(0) 



4r(l + m) 



(7) 



[Eq. (P12)]. 



L. Ciotti and B. Lanzoni: Stellar systems following the _R 1//m luminosity law. 



3 



3. The DF for anisotropic R 1 /™ models 

For any collisionless stationary system the DF / depends 
on the phase-space coordinates only through the isolating 
integrals of motion admitted by the underlying potential 
(Jeans Theorem, Chandrasekhar 1942), and moreover, if 
the system is also spherically symmetric in all its prop- 
erties, / depends only on the binding energy and on the 
angular momentum square modulus I? . Usually, the nega- 
tive value of the binding energy, the relative binding energy 
£, is used. For spherical models with / = f(£,L 2 ), the 
tangential components of the velocity dispersion tensor 
are identical, the only possible difference being between 
of and of = cr| = of /2, and the total velocity dispersion 
is er 2 (r) = of (r) + a t ( r )- 

In the OM formulation the radially anisotropic case 
is obtained assuming a DF depending on E and L 2 only 
through the variable Q defined as: 



T 2 
2ri 



(8) 



where r a is the so-called anisotropy radius, and where 
f(Q) = for Q < 0. Under this assumption the mod- 
els are characterized by radial anisotropy increasing with 
the galactic radius, and 



2o?(r) 



r 2 + r. 



2 ' 



(9) 



In the limit r a — > oo the velocity dispersion tensor is glob- 
ally isotropic. The simple relation between energy and an- 
gular momentum prescribed by Eq. (8) allows to express 
the DF as: 



f(Q) 



where 



1 



dg dip 



0(r) = 1 + -o )p(r 



(10) 



(11) 



(BT, p. 240). For ease of comparison we will use in the fol- 
lowing s a =jjRe, and f m (Q) = [G 3 T J i?;!]-i/2 x J^q^ 
where < Q = Q/ip m {0) < 1. 

3.1. Consistency 

The basic requirement for any physically admissible DF is 
its non-negativity over the phase-space accessible to the 
system, and we call consistent any model with a nowhere 
negative /. The DF obtained by Eq. (10) is not automat- 
ically consistent, because the integral inversion does not 
guarantee its positivity, and so for any model one has to 
check the consistency: if for some positive value of Q it re- 
sults / < 0, the adopted anisotropy radius is inconsistent 
with the assumed density profile. For example, in Paper 



I it was shown that globally isotropic R}l m models are 
consistent for all the explored values of m. It is therefore 
of interest to investigate here the consistency of the same 
family of models Jor various degrees of anisotropy. In Fig. 
1 we plot f m vs. Q in the case of global isotropy and strong 
anisotropy. A common characteristic is that the central di- 
vergence of the DFs, present in all globally isotropic mod- 
els, is unaffected by OM anisotropy, a behaviour similar 
to that analytically discussed in Ciotti (1996), and for 
which the same qualitative explanation holds. This diver- 
gence - as for any density profile with total finite mass - 
is not a problem: although the central phase-space den- 
sity of the models diverges, the corresponding mass does 
not. The requirement of consistency leads to define the 



Fig. 1. The DFs for globally isotropic (solid lines) and stable 
anisotropic (£ = 1.7, dotted lines) R 1 /™ models. 



critical anisotropy radius for consistency (s a ) c so that for 
Sa < (s a ) c a negative DF for some admissible value of Q 
is obtained. The existence of a critical anisotropy radius 
for the R}/ m models is easily understandable, remember- 
ing that a completely radial orbital distribution cannot 
be sustained by density profiles less divergent than 1/r 2 
(see, e.g., Ciotti & Pellegrini 1992 for an easy proof), and 
that the logarithmic slope for R 1 / 771 models is 1/m—l for 
r — ► 0. It is then interesting to show in detail the effect of 
a decreasing s a on the DF of R l ' m models. In particular, 
in Fig. 2 we show the modifications on the de Vaucouleurs 
DF. Note that approaching (s a ) c the DF becomes more 
and more depressed in the regions^corresponding to inter- 
mediate values of the parameter Q, i.e., fa becomes first 
negative outside the center. The dramatic effect of the 
anisotropy when s a is near its critical value is apparent in 
Fig. 2, where we have plotted fa also for a slightly higher 
value of s a . This behaviour is common to the whole family 
of the R 1 / 771 models, and seems to be more a consequence 
of the OM parameterization itself than a characteristic of 
some specific mass model.El 

In Fig. 3 (solid line) (s a ) c is plotted in the parameter 
space (m; s a ). Note how (s a ) c asymptotically decreases to- 
wards very small values for increasing m. The qualitative 
trend of (s a ) c is due to the behaviour of v m : for small m 
this results in a quite flat density distribution, and so only 
"high" values for s a are permitted; the opposite is true for 
high m models, for the density becomes more and more 



1 An identical behaviour it is found also for one and two 
component Hernquist models, discussed in Ciotti (1996). 
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Fig. 2. The modifications of the de Vaucouleurs DF, moving 
from the globally isotropic case (solid line), to the £ = 1.7 
anisotropy (dotted line), and to the critical anisotropy for con- 
sistency (dashed line). 

similar to a profile oc 1/r, and a stronger radial anisotropy 
is admitted. The flattening of the curve (s a ) c for high m 
is explained by the same argument. As can be seen in Fig. 



4), according to the previous discussion. So, assuming a 
fiducial critical value of £ for stability, a minimum value 
for the anisotropy radius (s a ){ is obtained, i.e., all models 
with s a < (s a )| arc unstable (see also point 4 in Sect. 6). 

In Fig. 3 (s a ){ is shown for two different £, and the 
corresponding values are given in Tab. 1. As it is intu- 
itive, for each m, (s a )j > (s a ) c : all models in the strip 
(s a ) c < s a < (s a ){ are consistent but unstable. As in the 
consistency analysis, also for stability an increase of m 
corresponds to a decrease in (s a )^: this is due to the fact 
that inside s a the orbital distribution is nearly isotropic, 
and with increasing m a higher fraction of the mass is con- 
tained in the central regions of the model, thus exerting 
a more efficient stabilizing influence on the system (Poly- 
achenko 1987). A fit of the minimum anisotropy radius 



Fig. 3. The minimum value of the anisotropy radius for the 
model consistency (solid line), and that for the model stability, 
for two different values of £. The lines are the interpolating 
functions given in the text. 

3, a good fit of (s a ) c as function of m is given by: 

(s a ) c ~ e-°- 93m (1.32 - 9.85 l(T 3 m 2 + 3.28 l(T 3 m 4 ), (12) 

while the exact values are given for integer m in Table 1. 

3.2. Stability 

A given density model is not useful for applications on 
data-modelling if unstable. In Paper I it was shown in 
a rigorous way that globally isotropic B}l m models are 
stable. Unfortunately for anisotropic models the same 
approach is not possible, and so only approximate re- 
sults can be obtained (unless one performs a much more 
complex linear stability analysis). Here, as in Carollo, 
de Zeeuw, & van der Marcl (1995, hereafter CZM), the 
stability of anisotropic models is investigated in a semi- 
quantitative way using the radial-orbit instability indica- 
tor £ = 2K T /K t (see, e.g., Fridman & Polyachcnko, 1984), 
where K t = 2n J pa 2 r 2 dr and K r = 2ir J pa 2 r 2 dr are the 
tangential and the radial kinetic energies, and have been 
numerically computed. This parameter is known to be a 
robust indicator, i.e., it is quite independent of the as- 
sumed density distribution profile, and when £>1.5 -j- 2 
the model is likely to be unstable. Note that for any glob- 
ally isotropic model £ = 1, because 2K r — K t , while in 
presence of radial anisotropy 2K r > K t , and so £ > 1. 
For the investigated models £ = £(s a ,ra), and for a fixed 
m it decreases towards unity for increasing s a (see Fig. 



Fig. 4. The value of the stability parameter £ for various m 
and for an increasing anisotropy radius. 

for stability, as a function of m, and having assumed as 
stability requirement that £ = 1.7, is 

(s a )| ~ e-°- 019m2 (1.31 - 3.17 10- 4 to 2 + 1.06 1(T 4 to 4 ).(13) 
4. Velocity Dispersion Profiles 

In this section we present the spatial and line-of-sight ve- 
locity dispersions profiles for the i? 1 /" 1 models with radial 
orbital anisotropy, and we compare them with the analo- 
gous globally isotropic cases, fully described in Paper I. In 
order to obtain the radial component tr r (r), we integrate 
the Jeans equation 

1 dp(r)a 2 (r) | 2 /3(r)<r r 2 (r) ^ rfy.(r) 
p(r) dr r dr 

with the natural boundary condition pa 2 — > for r — > 
oo. Having assumed OM anisotropy, the integral solution 
can be written explicitly, as shown by Binney & Mamon 
(1982), and after normalization cr 2 (r) = GTI R C x ct 2 (s). 
The tangential velocity dispersion is then obtained from 
Eq. (9). In Fig. 5, <r r /<7y and ag/ay are shown for some 
values of m, where a v = (2n/M) J °° pipr 2 dr, is the virial 
velocity dispersion. After normalization it results that a 
very good fit is given by 

a 2 ~ 4.7e- 1 - 82m . (15) 

A common feature of all the models is the charac- 
teristic central depression of the velocity dispersion pro- 
files: the explanation of this behaviour for isotropic models 
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Table 1. Critical values of s a for consistency and stability in the cases £ = 1.7 and £ = 2, for integer values of m 
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Fig. 5. The isotropic (solid lines), radial (dotted lines), and one-dimensional tangential (dashed lines) velocity dispersion profiles 
for various m and anisotropy radii. For each model the s a is assumed to be the minimum possible for stability (with £ = 1.7), 
and its value is the number printed in each panel. 



was qualitatively given by Binney (1980) for the de Vau- 
couleurs law, and analytically for all the R}l m models in 
Paper I. We note that with increasing m the maximum 
moves towards smaller and smaller radii (for the m = 10 
model the maximum is inside s = 10~ 2 ), and its value 
becomes correspondingly higher and higher. The same be- 
haviour is shown also in the strongly anisotropic models, 
both in er r and a t . It can be noted how the position and 
the value of the maximum arc not strongly affected by 
anisotropy. 

4-.1. Projected and Aperture Velocity Dispersions 

More important for observational purposes is the line- 
of-sight (or projected) velocity dispersion profile <7p(i?), 
obtained from cr r as: 



I(R)al(R) = 2 



n>2 

r z 



u(r)a 2 (r)rdr 



\fr 



R 2 



(16) 



(see, e.g., BT, p. 208). 

In Fig. 6 <7p/<7v corresponding to the same models 
described in Fig. 5 is shown. Note that, as a consequence 
of projection, the central depression is somewhat reduced 
but does not completely disappear. We give here a simple 
approximation of the position of the maximum and its 
height, as a function of m, in the case of global isotropy: 



— ] ~0.62e' 



0.07m 



°V/ Max 



R 



Max 



Rp. 



e -°' 13m2 (0.36- 



2.2710 



" 3 m 2 



(17) 



-7.5810~ 4 m 4 ).(18) 



These fitting formulae may be useful when correcting the 
observed values of the velocity dispersion, in simple appli- 
cations of the virial theorem. 

When observed through an aperture of finite size, the 
projected velocity dispersion profile is weighted on the 
brightness profile I(R). As in Ciotti et al. (1996), we ap- 
proximate this quantity calculating the aperture velocity 
dispersion, defined as 



LP(R)al(R) = 2ir / I{R')al{R')R' dR' , 



(19) 



where LP{R) is the projected luminosity inside R [Eq. 
(P3)]. In Fig. 7 we plot a a (R)/ay for the same models of 
Fig. 5. Note how, independently of the anisotropy radius, 
c a — > ay/V3 for R — ► oo: this result can be proved to be 
true for any assumed anisotropy (see, e.g., Ciotti 1994). 

4-2. Implications on the FP 

Looking at our results on the velocity dispersion profiles, 
we are tempted to discuss qualitatively their implications 
on the problem of the FP tilt and thickness. 

Two main considerations can be made. The first con- 
cerns the effect of radial anisotropy as a possible origin of 
the FP tilt, through a systematic increase of it with galaxy 
luminosity. From Fig. 7, where the maximum degree of 
(OM) anisotropy consistent with stability is considered, 
it is clear that the radial anisotropy cannot produce - in 
the assumption of structural homology (i.e., the same m 
for all galaxies) - the required variation of a factor of 3 in 
the observed squared velocity dispersion. This conclusion 
was already reached for different galaxy models in Ciotti 
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Fig. 6. The projected velocity dispersion for the same models shown in Fig. 5. The solid line refers to the globally isotropic 
models, and the dotted line to the anisotropic ones. 



Fig. 7. The aperture velocity dispersion for the same models shown in Fig. 5. The solid line refers to the globally isotropic 
models, and the dotted line to the anisotropic ones. 



et al. (1996). In any case - at least in principle - a possible 
observational test for the importance of anisotropy would 
be to construct the FP using cr a measured at large radii, 
and see whether its tilt is reduced. 

The second point concerns the problem of the very 
small thickness of the FP. From Fig. 7 one can see 
that for increasing m the maximum deviation between 
the isotropic and anisotropic velocity dispersions becomes 
smaller and smaller: for example, the percentage difference 
of their squared values for m = 2 is ~ 15%, for m = 4 it 
is already reduced to ~ 7%, and for m = 10 is less than 
6%. This implies that the anisotropy is not required to be 
fine-tuned with the galaxy luminosity in order to maintain 
the small observed FP scatter (< 12% in o\; see, Ciotti et 
al. 1996); on the contrary, all the admissible anisotropies 
can be present at each luminosity maintaining at the same 
time the FP thin. 



5. Velocity Profiles 

The velocity profile (VP) at a certain projected distance 
from the galaxy center is the distribution of the stars line- 
of-sight velocities at that point. It is strictly linked to the 
line profile in the observed spectrum, that is the convolu- 
tion of the stellar spectra with a certain VP. The shape of 
the VPs depends not only on the galactic potential, but 
also on the stars orbital distribution, a dynamical property 
not fully determined by the galactic potential itself. That 
is why the usual assumption of a gaussian shape for the 
VPs can generate a loss of information, and for this reason 
it has been suggested that the analysis of the deviations 
of VPs from gaussianity may give important insights on 
the dynamical structure of a galaxy (van der Marel 1994). 

In the case of OM anisotropy the VPs can be numer- 
ically recovered from the f(Q) using the integral expres- 



sion given by CZM. As usual, we expand the VPs on the 
Gauss-Hcrmite basis: 



VP(v) 



7e2< 



N 



'2 it a 



J2hjHj(v/a), 



(20) 



3=0 



(Gerhard 1993, van der Marel & Franx 1993), where the 
Hj are the Hermite polynomials as given in van der Marel 
& Franx. The even coefficients h 2 j measure symmetric de- 
viations from a pure gaussian, while the odd coefficients 
/i2j+i are identically zero for our models because the DF 
depends only on L 2 . We limit our investigation to the coef- 
ficient hi, because higher order coefficients are usually not 
available from spectroscopic observations. The uncertain- 
ties in the published data on hi are of the order of <~ 0.02 
(CZM). In practice, we fix ho = 1, and h 2 = in Eq. 
(20), thus requiring that the first term is the best-fitting 
gaussian, and we minimize the \ 2 using the Levenberg- 
Marquardt method (Numerical Recipes, p. 678) for the si- 
multaneous non linear fit of (7, a, hi). In Fig. 8 the radial 
behaviour of a and hi for various m and for isotropic (left 
panels) and anisotropic (£ = 1.7, right panels) models are 
shown. 

Isotropic R}l m models have surprisingly gaussian VPs, 
as the smallness of their hi indicates (Fig. 8a). Exception 
is made by the low-m models, which present strongly non- 
gaussian VPs at small radii. Note also how, for R>0.7R e 
the hi are completely indistinguishable for all m and indi- 
cate (slightly) flat-topped VPs (as generally a negative hi 
indicates). Moving inside, the differences between models 
become more and more important, but the global trend is 
towards VPs more peaked than the best-fitting gaussian, 
especially for low values of m. The radial behaviour of the 
corresponding dispersions a of the best-fitting gaussian is 
shown in Fig. 8b, where the similarity (apart a re-scaling) 
with the dispersions plotted in Fig. 6 is evident. 

For anisotropic models the a (Fig. 8d) are very similar 
to that of the isotropic case, while the situation is quite 
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Fig. 8. The values of the parameters /i 4 (upper panels) and a (normalized to (GT/ _R ) 1/ ' 2 , lower panels) characterizing the 
VPs of the R}^ m models, as a function of R, for m = 1 (circles), m = 2 (crosses), m = 4 (squares), and m = 10 (triangles). Left 
panels refer to the isotropic models, right panels to anisotropic (£ = 1.7) models. 



different for the /14 (Fig. 8c). Their values systematically 
decrease with radius in the inner regions, but rapidly in- 
crease at large radii, indicating significantly top-peaked 
VPs. This general trend, and the values as well, are nearly 
the same as those found by CZM for the 7 models, and are 
due to the orbital distribution in the outer part of models 
with OM radial anisotropy. 

Because of the differences in the trend of /14 between 
the isotropic and the anisotropic cases, and because of the 
growing evidence that i? 1 / m -law appropriately describe 
the surface brightness profiles of elliptical galaxies, we 
conclude that a detailed study of the VPs along the FP 
could be in principle a tool to study the effect of orbital 
anisotropy on its tilt and thickness. 

6. Conclusions 

The results of this work are the following: 

1. For OM anisotropic i? 1 /" 1 systems the consistency re- 
gion in the parameter space is explored, and the mini- 
mum anisotropy radius that can be assumed for given 
to in order to have a physical model is determined. 
We find that models with high to can sustain more 
anisotropy than models with low to. A fit of the criti- 
cal anisotropy radius for consistency is given. 

2. The maximum anisotropy tolerated by R}l m models 
in order to be stable against radial orbit instability is 
approximately derived, with the aid of the standard 
stability parameter £. As expected, in this case the 
limitation on s a is stronger than that required by con- 
sistency. Again, high-m models are more stable than 
low to models. A fit of the minimum anisotropy radius 
permitted for stability is given. 

3. The spatial, projected, and aperture velocity disper- 
sions are derived for various degrees of anisotropy. 
Their main characteristic is the fact that the off-center 
maximum - a feature already known and discussed for 
globally isotropic R 1 / 711 models - is still present, and 
not very much affected even by a strong anisotropy. 

4. The implications of this work for the problem of the tilt 
and thickness of the FP of elliptical galaxies are that 
orbital anisotropy cannot be at the origin of the tilt 
if galaxies are described by the R}l m law and char- 
acterized by structural homology. At the same time, 
the small thickness of the FP at fixed luminosity does 



not require any fine tuning between anisotropy and 
luminosity, due to the stability requirement. We note 
that the Saha (1991, 1992) works on stability further 
strengthen our conclusions: he found that radial or- 
bit instability can affect also models with £ smaller 
than the values suggested by Fridman & Polyachenko 
(1984) and here used, and so a still smaller amount of 
anisotropy would be permitted. 
5. The VPs are studied at various distances from the cen- 
ter for different anisotropy degrees and values of m. For 
globally isotropic models the VPs are very well approx- 
imated by a gaussian, except for very small radii and 
for low to, where detectable deviations from a pure 
gaussian are revealed. The lower order correction, pa- 
rameterized by the coefficient hi, shows that outside 
R e the VPs are flat-topped and essentially indistin- 
guishable for different values of to. On the contrary, 
the VPs of anisotropic models at R>R e are more cen- 
trally peaked than a gaussian, and the values of /14 
increase with increasing to. 
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